1-25  il-4 


6 


— ■  A.  na  "it 


DISTRIBUTION  STATEMENT  /I 

Approved  for  public  release! 
Distribution  Unlimited  > 


At>-  A  192. 

Journal  of  Chemical  Physics,  in  press 


ANOMALIES  IN  THE  HEAT-CAPACITY  SIGNATURES  OF  SUBMONOLAYER  ADSORBATES  WITH 


ATTRACTIVE  LATERAL  INTERACTIONS 


Young  Sik  Kim1 
Department  of  Chemistry- 
University  of  Rochester 
Rochester,  New  York  14627 


Franco  Battaglia*  and  Thomas  F.  George 
Departments  of  Physics  &  Astronomy  and  Chemistry 
239  Fronczak  Hall 

State  University  of  New  York  at  Buffalo 
Buffalo,  New  York  14260 


DTIQ 

^LECTEI 


M  19 
09 

D 


Abstract 

The  analytic  closed  form  of  the  heat -'capacity  signatures  previously  derived 
for  the  McQuistan-Hock  (MQH)  model  of  a  lattice  gas  is  applied  to  various 
adsorbed  systems  for  which  the  lateral  interaction  varies  from  a  few  meV  to  about 
300  meV.  It  is  shown  that  whenever  the  adsorption  system  can  be  described  by  a 
two-dimensional  gas  on  which  the  substrate  effects  are  less  important  than  the 
adatom-adatom  interactions,  the  computed  temperatures  at  which  the  heat-capacity 
signatures  display  their  maximum  are  in  excellent  agreement  with  the  experimental 
measurements . 
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I .  Introduction 

Recently  McQuistan  and  Hock^  have  developed  an  exact  solution  for  the 

distribution  function  of  q  indistinguishable  particles  on  a  2  x  N  lattice,  where 

2 

N  is  a  positive  integer.  In  a  previous  paper  (hereafter  called  paper  I),  we 
have  derived  an  analytic  closed  form  for  the  heat-capacity  signatures  from  the 
McQuistan- Hock  (MQH)  model  and  have  tested  the  model  by  comparing  our  computed 

temperatures  T  at  which  the  heat  capacity  displays  its  maximum  value  with 

»  c 

experimental  measurements  on  Ne ,  Ar  and  Xe  adsorbed  on  graphite,  and  with  some 
other  model  calculations.  In  spite  of  its  simplicity  and  its  being  an  almost 
one  -  dimens ional  model,  the  MQH  model  gave  for  these  systems  good  predictions  for 
the  temperatures  at  which  the  heat-capacity  signatures  display  their  anomaly. 
For  instance,  denoting  by  R  the  ratio  between  Tc  and  the  interaction  strength 
|V^|  among  two  nearest -neighbor  occupied  pairs,  the  one -dimens ional  Ising 

*1  r\ 

model,  the  MQH  modelz  and  the  average  experimental  measurements  for  Ne ,  Ar  and 
Xe  on  graphite  give  for  R  the  values  of  0.21,  0.43  and  0.46,  respectively.  The 
MQH  model  obviously  represents  a  logical  extension  from  a  one -dimensional  to  the 
more  realistic  full  two-dimensional  representation  of  the  physical  situation. 

One  might  wonder  about  the  justification  of  even  considering  a  2  x  N  model 
given  that  there  are  excellent  numerical  techniques  for  finding  with  well- 
controlled  accuracy  the  solution  to  the  full  two-dimensional  problem.  The 
justification  resides  in  the  very  fact  that  in  going  from  the  one -dimensional 
model  to  the  2  x  N  one,  there  is  an  abrupt  significant  improvement  in  the 
predictions  of  the  experiments  as  far  as  the  Tc's  are  concerned.  This  implies 
not  only  that  moving  to  the  N  x  N  model  should  not  improve  the  results 
significantly,  but  also  and  more  importantly,  that  most  of  the  physics  of  the 
liquid-vapor  equilibrium  in  two -dimens ions  is  already  in  the  2  x  N  quasi-one- 
dimensional  model.  Actually,  as  we  shall  see,  the  model  describes  only  a  limited 
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portion  of  the  full  phase  diagram  of  a  two-dimensional  system,  namely,  the  one 
relative  to  the  liquid-vapor  equilibrium.  Furthermore,  the  model  has  the 
advantage  of  being  exactly  solvable  and  containing  very  few  parameters,  namely 
^11’  V0  (c^e  interaction  strength  between  an  adparticle  and  its  own  site)  and  V00 
(the  interaction  strength  between  two  vacant  nearest-neighbor  sites) .  This 
latter  parameter  might  turn  out  to  be  quite  valuable  in  that  it  accounts  for 

distortion  of  the  substrate,  an  effect  that  has  been  invoked  by  various 

4-6  . 

authors  to  explain  some  details  of  the  experimental  heat-capacity  signatures. 

Due  to  the  remarkable  agreement  between  the  experimental  values  of  and  the 

ones  calculated  from  the  MQH  model  for  the  systems  we  have  made  our  test,  we  have 

been  encouraged  to  apply  the  model  to  a  wider  variety  of  systems  so  as  to  scan  a 

greater  number  of  experimental  data.  Our  choices  have  been  limited  only  by  the 

availability  of  accurate  estimates  of  the  interaction  strengths  and  accurate 

experimental  heat  -  capacity  signatures.  In  particular,  we  have  considered  Kr ,  CH^ 

and  O2  on  graphite,  and  Cu,  Ag,  Au,  Ni  and  Pd  on  tungsten.  In  the  next  section 

we  review  briefly  the  MQH  model  in  connection  with  the  computation  of  heat- 

capacitv  signatures.  In  the  third  and  last  section  we  present  and  discuss  our 

results . 

II .  Theory 

The  theory  of  the  2  x  N  lattice  has  been  extensively  described  in  detail  in 
Ref.  1  and  its  application  to  the  computation  of  heat-capacity  signatures  in  Ref. 
2  (paper  I).  In  this  section  we  give  only  the  relevant  results.  In  order  to 
study  the  behavior  of  the  heat  capacity  as  a  function  of  temperature,  keeping  the 
spreading  pressure  p  and  the  number  of  adparticles  q  as  constant,  it  is  most 
convenient  to  work  within  the  isothermal - isobaric  ensemble,  whose  partition 
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function  for  a  system  of  q  indistinguishable  particles  on  a  2  x  N  lattice  is 
given  by 


Mq,0,p)  -  tq  ^  fN  (q,0)  r,N  , 


-  p  y 

where  rj  *■  e  ,  c  —  e  and  ^  ^  is  the  temperature  times  the  Boltzmann 


constant  k  .  In  Eq .  (1), 
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where  x  —  e  and  y  =  e  The  number  of  unique  ways  q  indistinguishable 
particles  can  be  arranged  on  a  rectangular  2  x  N  lattice  to  form  n^  occupied 
nearest-neighbor  pairs  and  ngg  vacant  nearest-neighbor  pairs  is  given  by 
A(N , q , Uqq , n^^ ]  by  making  use  of  a  15-term  recursion  relation.  In  paper  I  we 
have  shown  that  A(q,/9,p)  can  be  written  as 


A  -  )  y(z.)  z.  q  ,  0  <  q  <  2N 

q  L  J  j 
J-i 


where  the  z ^ ' s  are  the  solutions  of  the  equation 
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From  A(q,/9,p),  the  heat  capacity  is  given  by 

C(q.^.p)  -  V^A1  '  (f')2]  ’  (10) 

where  the  prime  denotes  the  derivative  with  respect  to  /3;  the  value  p  -  (kgTc) 
at  which  C(q,/9,p)  exhibits  its  maximum  is  given  by  the  solution  of  the  equation 

2A(A'  ) 2  +  3/9AA' A"  -  20(A')3  -  0&2A"  -  2A2A"  -  0  (11) 

III .  Results  and  Discussion 

The  MQH  model  contains  very  few  parameters,  namely  Vg,  Vqq  and  V^;  Yet 
those  parameters  are  apparently  sufficient  to  determine  the  "gross"  features  of 
the  statistics  of  adsorption.  In  this  paper  we  are  interested  at  the  simulation 
of  the  critical  temperatures  as  determined  by  the  heat-capacity  signatures. 
Since  the  effect  of  Vqq  is  expected  to  be  small  and  since  there  is  no  reliable 
estimate  for  it,  we  have  chosen  Vqq  =  0,  which  is  equivalent  to  assuming  a 

"rigid"  lattice  in  which  there  is  no  interaction  between  two  nearest-neighbor 
vacant  sites.  This  is  the  same  choice  made  by  Hock  and  McQuistan  in  their 
computation  of  adsorption  isotherms.2  Moreover,  as  can  be  seen  from  Eqs .  (1), 

(3),  (4) ,  (9)  and  (10),  although  the  partition  function  and  the  Gibbs  free  energy 
do  depend  on  Vq,  the  heat  capacity,  C(q,£,p)  in  the  thermodynamic  limit  does  not. 
The  role  of  the  substrate  is  analogous  to  that  of  a  third  body  in  the  description 
of  the  interaction  between  two  particles  embedded  in  a  three-dimensional  medium, 
and  the  only  task  the  substrate  accomplishes  is  to  force  the  adsorbed  particles 
onto  a  plane.  The  all-important  parameter  here  is  V^,  i  ®-,  c^e  interaction 


strength  between  two  occupied  pairs,  whose  value  is,  of  course,  different  than 

the  free-space  value.  For  Kr  and  CH^  we  refer  to  the  results  obtained  by  Cole 

and  coworkers^’^  which  are,  to  our  knowledge,  the  most  reliable  data  available. 

Their  value  of  jv^l  -  170  K  for  Kr  is  in  remarkable  agreement  with  Putnam's10 

value  of  171  K.  For  CH^  we  take  the  value  of  |v^|  ”  ^77  K.  Another  widely- 

studied  system  is  two-dimensional  oxygen  adsorbed  on  graphite.  It  has  been 

investigated  experimentally11  with  use  of  heat  capacity  measurements,  and 

12 

calculations  have  been  performed  on  it  especially  by  Etters  and  coworkers. 

13 

These  authors  have  used  (and  so  do  we)  the  estimate  by  English  and  Venables  for 

|V11|,  i.e.,  54  K.  All  systems  mentioned  so  far,  together  with  the  ones  studied 

in  paper  I,  are  adsorbed  systems  with  a  lateral  interaction  of  the  order  of  10 

meV.  Recently,  Kolaczkiewicz  and  Bauer1^  have  studied  chemisorbed  layers  of 

metals  (Cu,  Ag,  Au,  Ni  and  Pd)  for  which  thermal  desorption  spectroscopy 

measurements  have  shown  that  their  desorption  energies  from  W[100]  surfaces 

increase  with  coverage,  implying  that  the  lateral  interactions  which  are  on  the 

order  of  300  meV  are  attractive.  We  shall  compare  our  results  for  this  range  of 

adatom-adatom  interactions  with  the  ones  of  Kolaczkiewicz  and  Bauer,  thereby 

showing  that  the  MQH  model  has  in  fact  a  wide  range  of  applicability. 

Critical  temperatures  for  the  mentioned  systems  have  been  determined  by 

15-21  ,  .  , 

various  techniques,  especially  from  heat  capacity  signatures,  hign- 

no  o  o  24 

resolution  synchrotron  X-ray  ’  and  low  energy  electron  diffraction  studies, 
thermal  desorption  spectroscopy , ^ ^  and  from  the  temperature  dependence  of  the 
work  function  of  adsorbate -covered  surfaces.  As  shown  in  the  previous  section, 
we  are  able  to  extract  from  the  MQH  series  expansion  for  the  canonical  partition 
function,  (Eq.  (21)),  a  closed  analytic  form  for  the  isothermal- isobaric 
partition  function  (Eq.  (3))  from  which  one  can  compute  heat-capacity  signatures, 
(Eq.  (10))  whose  maximum  is  exhibited  at  the  temperature  T£  that  solves  Eq.  (11). 


In  the  MQH  model  the  dependence  of  Tc  on  is  essentially  linear  with  a  slope 

that  depends  smoothly  on  the  coverage  6  (in  particular,  in  the  range  0.2  <  9  < 

0.7,  the  slope  varies  between  0.37  and  0.44).  In  Figs.  1  and  2  we  display  Tc  as 

a  function  of  |V.^|  for  two  ranges  of  |V^|,  namely  0  -  300  K  (Fig.  1)  and  2000- 

4000  K  (Fig.  2).  To  compare  with  data  available  from  literature,  we  have  chosen 

9  -  0.5  for  the  systems  on  graphite  and  9-0.3  for  the  ones  on  tungsten.  With 

the  exception  of  Kr  and  Xe,  the  agreement  between  the  prediction  from  the  MQH  and 

the  experimental  values  for  the  systems  shown  in  Fig.  1  is  very  good.  According 
28 

to  Dash,  one  can  distinguish  between  three  temperature -dependent  types  of  film 
growing:  either  three-dimensional  bulk  develops  asymptotically  from  the  two- 
dimensional  system,  or  there  is  no  film  but  cluster  growth,  or,  after  the  first 
two  or  three  monolayers,  cluster  growth  begins.  Clearly,  systems  for  which  the 
film  grows  monolayer  after  monolayer  are  suitable  to  be  well  simulated  by  a  two- 
dimensional  lattice -gas  model  such  as  the  MQH  one. 

In  order  to  understand  why  the  model  gives  such  excellent  predictions  for 
Ne,  C>2 ,  Ar  and  CH^  on  graphite,  but  not  for  Kr  and  Xe ,  let  us  first  recall  that 
adsorbed  systems  whose  two-dimensional  space  symmetry  results  from  the  two- 
dimensional  space  symmetry  of  the  substrate  surface  by  adding  or  subtracting 
symmetry  elements  are  named  commensurate,  and  that  adsorbed  structures  whose 
symmetry  is  not  related  to  that  of  the  substrate  are  named  incommensurate.  In 
many  cases  the  periodic  substrate  of  the  hexagonal  basal  surface  of  graphite  is 
covered  by  a  /3  structure,  in  which  the  centers  of  non- adjacent  hexagons  are 
preferred  adsorption  sites  and  form  a  triangular  lattice.  Whether  an  adsorbate 
has  commensurate  phases  depends  on  its  incompatibility  i,  defined  as  i  -  (a-d)  x 
100/d,  where  a  is  the  lattice  constant  of  the  /3  structure  (4.263  A  for 
graphite^®)  and  d  is  the  lattice  constant  of  the  [111] -plane  of  the  adsorbate 
species  in  three-dimensional  bulk  (Ne,  Ar,  Kr,  Xe,  CH^  and  O2  condense  into  the 


cubic  system) .  For  incompatibilities  i  positive  and  less  than  a  substrate- 


dependent  value  i  ,  commensurate  structures  are  energetically  preferred.  For 
graphite  i  -  5 ,  and  the  values  of  i  for  Ar,  Kr,  Xe  and  CH^  on  graphite  are  14.6, 
5.4,  -2.8  and  1.6,  respectively . ^ ^  Adsorbates  in  a  commensurate  phase  cannot 
be  properly  described  by  a  lattice -gas  model  such  as  the  MQH  one,  in  which  the 
only  role  of  the  substrate  is  to  constrain  the  adsorbate  particles  on  a  plane. 
In  an  incommensurate  phase  the  adatom- adatom  interaction  dominates  on  the 
substrate  influence,  and  the  MQH  model  is  expected  to  be  well  suited  to  simulate 
the  heat-capacity  anomalies  corresponding  to  two-dimensional  incommensurate 
phases . 

From  the  above  values  of  the  incompatibilities  and  from  the  fact  that  there 

is  no  evidence  of  a  commensurate  phase  for  (for  which  the  incompatibility 

value  can  be  safely  estimated  to  be  i  -  10) ,  it  is  now  clear  why  the  MQH  model 

gives  an  excellent  prediction  of  the  heat-capacity  anomalies  for  0£  and  Ar,  and  a 

poor  prediction  for  Kr.  (Kr  displays  a  commensurate- incommensurate  phase 

transition  '  but,  in  the  range  70-130  K,  it  occurs  at  coverages  greater  than 

3  A-  35 

unity;  Ar  was  found  incommensurate  at  all  conditions  studied.  ’  )  As  already 

mentioned,  CH^  has  a  low  value  of  incompatibility  which  might  invalidate  the 
above  reasoning.  This  is  not  the  case  because  the  phase  diagram  of  two- 
dimensional  methane  adsorbed  on  graphite^  shows  a  temperature- triggered 
commensurate- incommensurate  transition  at  T  -  47  K  and  coverages  below  0.8 

*>  r 

monolayers.  Therefore,  the  anomaly  at  T  -  75  K  occurs  in  a  region  where  the 
adsorbate-adsorbate  interactions  play  a  major  role. 

This  is  also  the  case  of  Ne,  as  can  be  seen  from  the  phase  diagram  proposed 
by  Huff  and  Dash^  for  submonolayers  in  the  temperature  range  1-20  K:  melting 

occurs  at  13.5  K,  so  that  the  anomaly  at  16  K  cannot  involve  any  solid 

commensurate  phase.  (The  barely  stable  ( Jl  x  /7)R19°  structure  exists  only  at 


10 


temperatures  below  4  K.  )  Xe  has  a  negative  incompatibility  and  is  therefore 
expected  to  be  highly  mobile,  because  when  the  [111] -plane  is  isolated  from  the 
bulk,  its  lattice  constant  (that  in  bulk  is  already  larger  than  the  lattice 
constant  of  the  structure)  becomes  even  larger  and  the  incommensurate  phase  is 
generally  more  stable.  However,  xenon  has  already  shown  a  peculiar  behavior  with 
its  unusually  high  value  of  the  ratio  ^2c^3c  i-n  spite  of  its  negative 
incompatibility.^^  This  has  suggested^  that  some  different  mechanism  of 


localization  is  causing  the  formation  of  a  structure  which  is  in  registry  with 
the  solid  substrate  and  with  a  positive  value  (11.6)  of  the  incompatibility, 
whatever  the  case  might  be,  we  point  out  that  although  the  quantitative 
predictions  from  the  MQH  model  for  adsorbates  whose  structure  is  somehow 
reminiscent  of  the  substrate  structure  are  not  as  excellent  as  for  adsorbates 
that  are,  without  doubt,  incommensurate  to  the  substrate,  the  relative 
qualitative  behavior  is  in  agreement  with  the  experimental  findings,  and  the 
quantitative  differences  are  not  so  serious  in  spite  of  the  simplicity  of  the 
model.  This  is  clearly  shown  for  Kr  on  graphite  and  for  the  systems  to  be 
considered  next. 

In  Fig.  2  we  compare  our  results  from  the  MQH  model  in  a  range  of  much 
stronger  (200-300  meV)  lateral  interaction  energies  with  systems  for  which 
accurate  estimates  of  both  adatom-adatom  potentials  and  critical  temperatures  are 
available.  Kolaczkiewicz  and  Bauer^  have  filled  a  gap  in  the  literature  by 
studying  chemisorbed  layers  with  attractive  lateral  interactions.  They  have 
constructed  from  work  function  measurements  the  coexistence  line  of  Cu,  Ag,  Au, 
Ni  and  Pd  adsorbed  on  W[110],  and  have  shown  that  at  submonolayer  coverages  these 
systems  can  be  described  by  a  law  of  corresponding  states  in  the  form  of  a  two- 
dimensional  van  der  Waals  equation  (with  a  critical  coverage  at  6  =■  0.3).  From 
their  adatom-adatom  Lennard- Jones  potential  parameters,  we  estimate,  taking  into 


•  • 
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account  the  harmonic  zero-point  energy  correction  2  J t / n)  ,  the  interaction 

strength  (a  is  the  hard-core  diameter  of  the  adatom,  e  is  the  well  depth  and 

p  is  the  reduced  mass).  This  is  displayed  in  Table  1  with  the  temperatures  Tc 

calculated  from  the  MQH  model  and  estimated  from  a  van  der  Waals  fitting  of  the 

experimental  data.  Again  we  can  see  that  the  qualitative  behavior  of  Tc  as 

computed  from  the  MQH  model  is  in  agreement  with  available  data  from  the 

literature,  even  in  the  region  of  strong  lateral  interactions.  The  van  der  Waals 

fitting  performed  by  Kolaczkiewicz  and  Bauer  on  their  experimental  data  gives 

interaction  strengths  that,  when  used  in  the  MQH  model,  give  back  the 

experimental  T^  trend.  Quantitatively,  the  results  are  not  as  satisfactory  as 

for  Ne ,  Ar,  0^  and  CH^  on  graphite.  We  cannot  ascribe  the  discrepancies  to  a  bad 

estimate  of  the  adatom-adatom  interaction  strength  that,  for  the  systems  on 

tungsten,  has  not  been  determined  as  accurately  as  for  the  systems  on  graphite. 

On  the  contrary,  we  expect  that  as  given  in  Table  I  is  an  upper  bound  to  the 

real  value,  thereby  giving  a  less  satisfactory  result  for  T£ .  Again,  the  reason 

for  this  has  to  be  found  in  the  structural  properties  of  these  metal  layers  on 

,  .  ,  26  ... 

tungsten  surfaces:  the  lateral  periodicity  of  the  adsorbate  is  very  close  (it 
not  identical)  to  that  of  the  substrate  whose  effects  on  the  Tc's  are  noC  taken 
into  account  by  the  MQH  model. 

In  Fig.  3  we  show  the  two-dimensional  phase  diagram  which  can  be  obtained 

within  the  limits  of  the  MQH  model.  The  displayed  curve  divides  the  plane  into 

two  regions  that  are  very  much  reminiscent  of  the  liquid-vapor  equilibrium  of  the 

bulk  phase.  Compared  with  available  experimental  data  of  chemisorbed  metals  on 

tungsten,  we  find  that:  (1)  the  critical  coverage  (©c  =  0.8)  in  the  MQH  model  is 

much  larger  than  the  experimental  results  (©c  =*  0.3),  and  (2)  the  ratios  -Tc/V^ 

2 

for  chemisorbed  metals  are  larger  than  those  from  the  MQH  model.  Since  this 
model  contains  only  the  interaction  strength  between  two  occupied  pairs  (V^) • 


-  .  «■  .  •  ,  •  .  -  _  *  .  •  «  "  _  -  .  A  ^  _%<  •  *>  w  w  *•  W V  m  '  ta^ 
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does  not  take  into  account  the  adparticle  size  as  well  as  substrate  effects.  If 
we  assume  the  simple  Lennard- Jones  pairwise  interaction  potential  between 
adparticles  and  fix  the  monolayer  density,  then  ©c  is  inversely  proportional  i.o 
the  hard-core  diameter  a  of  the  adatom.  Thus,  the  MQH  model,  which  neglects  the 
size  of  the  adatom,  will  give  large  values  for  the  critical  coverage.  Also,  the 
adatoms  are  bound  more  strongly  to  the  substrate  atoms  than  to  themselves,  so 
that  fewer  electrons  are  available  for  lateral  bonding.  The  strong  adatom- 
substrate  interaction  reduces  V^,  which  in  turn  increases  the  ratio 
should  further  mention  that  at  very  low  and  at  high  coverages  the  model  might 
fail  because  of  the  substrate  inhomogeneity  and  because  of  multilayer  effects, 
respectively . 

In  conclusion,  the  MQH  model  is  very  appealing  for  describing  heat-capacity 

anomalies  of  adsorbate  -  substrate  systems,  provided  (i)  no  solid  commensurate 

phase  is  involved  in  the  region  around  Tc>  (ii)  che  pairwise  interaction 

among  nearest-neighbor  adparticles  is  the  parameter  that  plays  the  major  role  in 

transition  process,  with  only  lower-order  effects  due  to  the  substrate,  and  (iii) 

the  transition  is  not  of  first  order.  These  conditions  are  not  very  restrictive, 

and  moreover,  they  are  not  independent  because  very  often  the  first  one  implies 

the  second,  as  has  been  already  discussed,  and  in  many  instances  also  the  third 

one,  since  several  transitions  from  or  into  commensurate  phases  have  been 

38 

predicted  to  be  of  first  order. 
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Table  I 


o 


Compuced  (T^)  and  experimental  (^£Xp)  critical  temperatures  of  Ne ,  Ar,  Kr ,  Xe ,  0^ 
and  CH^  on  graphite  (.coverage  0  -  0.5)  and  Cu,  Ag,  Au,  Ni  and  Pd  on  W[110] 
(coverage  8  =0.3).  are  the  adatom-adatom  interaction  strengths.  All  data 
are  in  Kelvin. 


-va) 

11 

T 

c 

& 

Ne 

35 

15 

16 

Ar 

120 

51 

50 

Kr 

170 

71 

85 

Xe 

236 

102 

118 

°2 

54 

23 

25 

CH4 

177 

75 

75 

Cu 

2850 

1120 

1170 

Ag 

2360 

927 

980 

Au 

2910 

1143 

1130 

Ni 

3430 

1348 

1400 

Pd 

2900 

1140 

1170 

a)  Taken  from  Refs.  8  and  9  for  Ne ,  Ar,  Kr,  Xe  and  CH^,  from  Refs.  12  and  13  for 

O2 ,  and  from  Ref.  14  for  Cu,  Ag,  Au,  Ni  and  Pd. 

b)  Taken  from  Ref.  8  for  Ne ,  Xe  and  CH^ ,  from  Ref.  21  for  Ar ,  from  Refs.  18-20 

for  Kr,  from  Ref.  16  for  O2,  from  Ref.  14  for  Cu,  Ag,  Ni  and  Pd,  and  from 

Ref.  27  for  Au. 


s  „>  .V,1-  A  «■  .  1A  . V.  -  A  AAAaAA  A  A.A  .V  A  A  A  A  A  A  /•  A  Av  A  A  AA  A  AA  -  ’jAy>  A  A 


1.  The  solid  line  shows  critical  temperatures  T£  from  the  MQH  model  as  a 


function  of  the  adatom-adatom  interaction  strength  V^,  in  the  range  0-30 
meV.  Circles  represent  available  data  for  some  adparticles  on  graphite. 


2.  The  solid  line  shows  critical  temperature  from  the  MQH  model  as  a  function 


of  the  adatom-adatom  strength  V.. 1  in  the  range  150-350  meV .  Circles 


represent  results  from  Ref.  14. 


3.  Phase  diagram  from  the  MQH  model.  Solid  line:  model  phase  boundary;  -fc:Cu, 


A:Ni,  • :Ag,  Q:Pd,  +:Au  on  W(110)  from  experimental  data  (Ref.  14). 
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